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ABSTRACT 

We simulate the collapse of isolated dwarf galaxies using SPH + N-Body simulations including a 
physically motivated description of the effects of supernova feedback. As the gas collapses and stars 
form, the supernova feedback disrupts enough gas to temporarily quench star formation. The gas 
flows outward into a hot halo, where it cools until star formation can continue once more and the 
cycle repeats. The star formation histories of isolated Local Group dwarf galaxies exhibit similar 
episodic bursts of star formation. We examine the mass dependence of the stellar velocity dispersions 
and find that they are no less than half the velocity of the halos measured at the virial radius. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: dwarf — methods: N-Body 
simulations 



1. INTRODUCTION 

Dwarf galaxies have shallow potential wells. They are 
therefore subject to many baryonic processes that are 
negligible in high mass galaxies. For example, dwarf 
galaxies are prone to gas loss through supernova feed- 
back and/or stripping. They also are unable to efficiently 
accrete gas whose temperature is higher than the virial 
temperature of the galaxy. As a result, the baryonic con- 
tent of dwarf galaxies is expected to be lower than the 
cosmic mean, and should have significant scatter. 

Given their susceptibility to galactic "weather" , the 
cool baryonic content of dwarf galaxies should fluctuate, 
not just from galaxy to galaxy, but within an individ- 
ual galaxy as a function of time. The expected temporal 
variations in the amount of cool gas should be accom- 
panied by variations in the star formation rate (SFR). 
These variations in the SFR can be seen within the Local 
Group, where dwarf galaxies are sufficiently well resolved 
that their past SFRs can be directly inferred from their 
color-magnitude diagrams. 

As expected, past star formation in dwarf galaxies ap- 
pears to be complex, with many sys tems e xhibitin g mul- 
tiple episodes of star forma tion (e.g. i Smecker-Hane et al 
TMe) [ Dohm-P almer et al.||20 02 Riz zi et al.||i^004| Itjki 



man|2005| jMcConnachie et al.|2006i ) . The temporal sep 
aration of bursts of star formation ranges from tens or 
hundreds of megay eaxs in dwarf irregulars like Sextans 
A, GR8 fc Phoe nix jDohm -Pal mer et al.J1998[[2002l[Dor 
phin et al.|2003 Young et al. 200 7|) , to gigayears m d warf 
spheroidals like Carina (Hurlcy-KeUeretaT 1998). In 
some systems these fluctuations have large amplitudes 
(e.g. Carina), and in others, the amplitudes are of the 
order of 50% of the mean star formation rate. However, 
the true instantaneous fluctuation amplitudes are pos- 
sibly larger if the fluctuation is shorter than the width 
of the time bins used for analysis. The temporal resolu- 
tion of star formation history recovery methods worsens 
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with increasing lookback time, so it is possible that short 
timescale fluctuations are present at early times as well. 

These fluctuating star formation rates have long been 
a puzzle. Many of the galaxies observed to have episodic 
star formation also have large intermediate age popula- 
tions, showing that they initially formed stars, then lay 
largely dormant, before reviving several gigay ears later 
(e.g. C arina, DDO210, LeoA, Leol, IC1613; ISmecker- 
Hane e t al. 1996: Hurley- Keller et al. 1998; McConnachle 
et al.||2006i ,Tolstoy et a l. 1998; Gall art et al.,|1999, ^er- 



nandez et ai.|2000 Skillman et al. 2003|). Others have no 



sign of current star formation, but were actively forming 
stars wi thin the past 100 Myrs ( e.g. DDO210, Dp 06, 
Fornax; [McConnachie et al.||2006| [Weisz et al.]|2007| . 

How is it possible that star formation in dwarf galax- 
ies can stop and start over a range of timescales? Some 
of these systems are relatively isolated (DDO210, LeoII, 
IC1613), ruling against interactions as a driving force. 
Some authors have speculated that reionization may have 
held off significant SF until intermediate ages 



je.g. 



Bul- 



locket al. 2000: Skillma n et al.|2003||G"nedin fc Kravtsov 



2006| Read et al.||2006|), but such a mechanism makes 



it hard to explain eit her the significant pop ulations of 



i tjj^m^ truly old stars ^ tetson et al.|1998| |Buonanno et al.| 
, 1 ' IGrebel &: Harbeck„ 2003| ) or the shorter timescale n 



1998 



luctu- 



ations seen at the present day. However, there is growing 
evidence from dwarf spheroidal age-metallicity relations 
that gas accretion and outflows may have played a role, 
leading to G- and K-dw arf problems compara ble to those 



seen in the Milky Way ( jKoch et al.||2006b|"a l 

In this paper, we explore the role that self-regulation of 
feedback can play in driving the observed fluctuations in 
the star formation rate, particularly on short timescales. 
To do so, we use high resolution numerical simulations of 
isolated galaxies, including self-consistent star formation 
and SN feedba ck recipes calibrated to match Milky Way 
sized galaxies ('Stinson et al. 2006). We maximize reso- 
lution and minimize the effect of interactions by study- 
ing the evolution of uniform, collapsing spheres. These 
simulations are fully self-consistent, and thus differ from 
previous efforts to model supernovae explosions in dwarf 
galaxies, most of which involve implanting an energetic 
source into the galaxy and tracking the resulting winds 
( |Dekel fc Silk|1986| [BiTrkert fc Ruiz-Lapuente|1997l [Mac 
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Low & Ferrara|1999| Scannapieco et al.|2000[ 2001 |Rec- 


chi et al.ll2001 


|Efstathiou||2000 |Fragile et al.||2003| Fu- 


jita et al.l[2003 


IPelupessy et al.|2004||Hensler et al.|2004| 


Marcolini et a 


.|2006|). 



logical context sacrifices some details of the mass accre- 
tion history, it also allows our models to have suSiciently 
high resolution to capture the more relevant physics of 
the star formation and feedback processes. The resulting 
smoothly collapsing halos are sufficiently simple so that 
we can isolate features of the galaxies' behavior that are 
due solely to the details of star formation and feedback. 

With these simplified but self-consistent models, we 
show that star formation in low mass galaxies can un- 
dergo a natural "breathing" mode, where episodes of star 
formation trigger gas heating that temporarily drives gas 
out of the cool phase and into a hot halo. Subsequent 
gas cooling and gas accretion then allows gas to settle 
back into the halo, and star formation begins again. This 
episodic mode of star formation is superimposed on the 
net infall of gas into the galaxy, leading to galaxies with 
both episodic star formation and significant intermediate 
age populations, in line with observations. 

In f|2] we describe the smoothed particle hydrodynam- 
ics (SPH) code, galaxy models, star formation, and su- 
pernova recipes used for the simulations. S|3] shows the 
sequence of events that lead to the episodic star forma- 
tion we see in our simulations. ^4] compares the star 
formation in the simulations with that observed in Local 
Group dwarfs. 

2. METHODS 

To simulate the star formation history in dwarf galax- 
ies, we creat ed five exam ple galaxies in hydrostatic equi- 
librium with |Navarro et al | ( |1997| (hereafter NFW) den- 
sity profiles. They vary in total mass from IO^Mq to 
1O^"M0 and have velocities at the virial radius (u2oo) of 
10, 15, 20, 25, and 30 km s^^ comparable to Local Group 
Dwarfs. We do not consider halos smaller than IO^Mq, 
which do not have enough mass to retain gas reionized 
by the UV radiation from the first stars and AGN. 

Although examining galaxies outside of the full cosmo- 
logical context sacrifices some details of the mass accre- 
tion history, it also allows our models to have sufficiently 
high resolution to capture the more relevant physics of 
the star formation and feedback processes. The resulting 
smoothly collapsing halos are sufficiently simple so that 
we can isolate features of the galaxies' behavior that are 
due solely to the details of star formation and feedback. 

2.1. Models 

We created combined gas an d dark matter vi r ialized 



halos using the specifications of Kaufmann et al. (20071. 
The halos start with 10^ gas particles and 10^ dark mat- 
ter particles, a resolution at which numerical losses of 
angula r momentum become small in [Kaufmann et aT| 
( 2007 1 . Equilibrium NFW dark matter halos with a con- 
8 were created following [Kazantzidis et aL\ 



centration c 



(20041, where c 



r/rs 



is the virial radius and r. 



is the scale radius. The dark matter extends beyond the 
virial radius to keep the halos in equilibrium. Dark mat- 
ter inside the virial radius is at slightly higher resolution 
than outside such that 90% of the dark matter particles 
are inside the virial radius, but only 72% of the mass is. 



Gas was added to the halo following the same NFW 
density profile as the dark matter, but comprising only 
10% of the total halo mass. Each gas particle was as- 
signed an initial temperature to keep it in hydrostatic 
equilibrium with the halo before the gas begins radia- 
tively cooling. All gas particles started with equal mass. 
The force resolution is based on a fixed fraction of 
the virial radius (r^ir), so that the softening is set to 
r„i^/2000. 

Dark matter particles are left with their random equi- 
librium velocities while the gas is spun with a uniform 

circular velocity so that A = is^-AEl ^ _ q Q3g^ where jgas 

is the average specific angular momentum of the gas, E 
and A/ are the total energy and mass of the halo. No cos- 
mological Hubble fiow is included in the initial velocities. 

2.2. Numerics 

The simulations in this p aper were run using t he paral- 
lel SPH code GASOLINE (Wads ley et al.|2004| . GASO- 
LINE solves the equations of hydrodyna mics, and in- 
cludes the radiative cooling described in |Katz et al.| 
(19961. The cooling assumes ionisation equilibrium, an 
ideal gas with primordial composition, and solves for the 
abundances of each ion species. Th e scheme uses th e 
coUisional ionization rates reported in Abel et al. (|1997 1, 
the radiative recombinat ion rates from ljlack"( [T98l[ ) and 
Verner fc F erlancT ('1 996[ ), bremsstrahlung, ana line cool- 
ing tromjCen ( 1992 ) ."The energy integration uses a semi- 
implicit stitt' integrator independently for each particle 
with the compressive heating and density (i.e. terms de- 
pendent on other particles) assumed to be constant over 
the timestep. 

Stars are formed and feedback is calculated every 1 
Myr. The star formation and feedback recipes are the 
"blast wave model" described in detail in Stinson et al.l 
(20061, and they are summarized as follows. Gas parti- 
cles niust be dense (n„iin = O.lcm"^) and cool (T„iax = 
15,000 K) to form stars. A subset of the particles that 
pass these criteria are randomly selected to form stars 
based on the commonly used star formation equation. 



dt 



td 



(1) 



yn 



where M^, is mass of stars created, c* is a constant star 
formation efficiency factor, Mgas is the mass of gas cre- 
ating the star and tdyn is the gas dynamical time. This 
equation is integrated on 1 Myr star formation time 
scales in all of the simulations presented here. T he con- 
stant param eter, c* is tuned to 0.05 to match the |Kenni- 



cutt 



( 1998 1 Schmidt Law in the simulated Isolated Model 
MilEyWay used in jStinson et al.| (|2006| . 

At the resolution of these simulations, star particles 
represent a large group of stars (300 Mq). Thus, each 
particle has its stars partitioned into mass bins based 
on t he initial mass function presented in |Kroupa et al.| 
'1993 ) . These masses are correlated to stellar lifetimes as 



(pg; 

desci 



escribed in Raiteri et al. ( 1996 1. Stars larger than 8 Mq 
explode as supernovae during tne timestep that overlaps 
their stellar lifetime after their birth time. Similarly, 
Type la supernovae are included for stars smaller than 
8 AIq. The explosion of all supernovae is treated using 
the analytic mod el for blastwaves presen ted in [McKee 



& Ostriker ( 1977 1 as described in detail in Stinson et al 
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220 Myr 



290 Myr 



345 Myr 



515 Myr 



670 Myr 



750 Myr 



830 Myr 



935 Myr 



1055 Myr 



( 2006 1 . While the blast radius is calculated using the full 
energy output of the supernova, only half of that energy 
is transferred to the surrounding ISM, Esn = 5 x 10^*^. 
The other half of the supernova energy is radiated away. 

3. RESULTS 

We ran simulations of the five isolated, low mass, virial- 
ized galaxy halos for 10 Gyr (approximately 30 dynam- 
ical times). In the lowest mass galaxies, we observed 
that the star formation varie s episodically. Oth er groups 



have noted similar behavior. Dong et al. ( 2003 ) produced 



episodic star formatio n in one dimensional simulations. 



Pelupessy et al. ([2004| starts SPH simulations at an ad- 
vanced stage wnere tnere is a significant population of 
stars, and all the gas is in the disk. Using SN feedback 
similar to what is employed here, Pelupessy et al._ (2004| 
still note episodic star formation in model galaxies with 
similar mass to our 15 km s~^ halo. 

The observed fiuctuations are driven by supernova 
feedback. As gas accretes and stars form, the result- 
ing feedback is sufficient to resist the collapse of gas and 
drive it outwards. The subsequent drop in gas density 
shuts off the burst of star formation. However, super- 
nova remnants cool and gas can once again accrete and 
another wave of stars can form. We now discuss this 
behavior in more detail. 

Figure T] shows the sequence of events illustrating the 
episodic behavior in a single galaxy. First, the infall of 
cooling gas triggers a central star burst (shown in the left 
panels). The subsequent supernova feedback contains 
enough energy to eject gas from the shallow potential 
well of the small galaxy halo (center panels) . As the gas 
expands adiabatically, the density of the gas in the center 
drops below the star formation threshold (n = 0.1 cm~^) 
(right panels). The cessation of star formation activity 
then allows gas to cool, contract and again become dense 
enough to form stars. The process then repeats in an- 
other cycle of star formation, expansion and cooling. In 
this way, an episodic star formation history is created in 
an isolated galaxy without triggering by external forces. 
Timescales and amplitudes for this episodic behavior are 
given in Table [l] as a function of galaxy mass. 

Figure [2] shows the star formation histories of the four 
most massive halos. The lowest mass halo is excluded be- 
cause its gas cannot cool efficiently and forms few stars. 
The star formation clearly ffuctuates in the 15, 20, and 
the 25 km s~^ halos. In these isolated, idealized galaxies, 
the ffuctuations are strictly periodic, though they would 
be unlikely to be so in a full cosmological context. The 
periods of their oscillation are 292, 342, and 387 Myr, 
respectively. These times are characteristic of the dy- 
namical time of the central region of constant density 
gas that forms after the starburst (250 Myr for the 15 
km s~^ halo). The dense core of the galaxies can cool in 
less than 100 Myr, which is much shorter than the dy- 
namical time. Thus cooling does not significantly impede 
the collapse of the galaxy and formation of stars and thus 
the period is instead dominated by the dynamical time 
for the gas to contract. 

For larger galaxies, star formation never completely 
ceases due to the reduced effectiveness of feedback in 
deeper potential wells. Some episodic behavior is seen, 
superimposed above a constant level and with a slightly 
longer period. This constant star formation exerts an 



Fig. 1. — The sequence of events that produces episodic star 
formation in the 15 km s~^ run. Purple dots represent the stars. A 
star burst in a dense gaseous region (colored bright red) causes the 
outflow of gas that prevents more stars from forming for a period 
of time until the gas is able to cool once more. The cycle repeats 
itself as there is a new star burst that drives gas and ends star 
formation. 
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Fig. 2. — The star formation history of four dwarf galaxies. The 
star formation histories vary as the mass of the galaxy increases, 
so that there are fewer, more widely separated fluctuations in more 
massive galaxy potentials. 

outward pressure on the gas, slowing the rate of infall and 
cooling of the gas, and increasing the separation between 
bursts. 

In addition to the slight change in period, the ampli- 
tude of SFR variations decreases with mass. At higher 
masses, the supernova energy is a smaller fraction of the 
galaxy's gravitational potential. 

4. COMPARISON WITH OBSERVATIONS 

To facilitate comparison with observations, we have 
translated the SFHs shown in Figure [2] into the form 
most commonly derived from observed color magnitude 
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Fig. 3. — The star formation rate as a function of stellar age 
where the age of the stars has been averaged over log bins as a 
rough comparison to observations. The newest stars on the left 
have the smallest bins to reflect observations where recent stellar 
ages are more easily resolved than old ages. 



TABLE 1 
Simulation data 



Mtot 


V200 


(7„ 






SFR 


(M©) 


(km 


(km 


(M©) 


(Myr) 


Amplitude 




s-i) 


s-l) 








3.18 X 10* 


10 


4.5 


1.23 X 10* 


292 


4.17 


109 


15 


7.4 


2.34 X 10'^ 


292 


2.98 


2.5 X 10^ 


20 


11.4 


2.15 X 10^ 


342 


2.25 


5 X lO'' 


25 


15.1 


7.86 X 10^ 


387 


0.73 


8.6 X 10^ 


30 


20.1 


2.20 X 10** 




0.42 



It is the most realistic converter from simulation to obser- 
vational data available today. The absolute magnitudes 
of the model galaxies agree well with the real low mass 
galaxies (e.g. 15 km hal o has My=-1 0.7, comparable 
to Phoenix with My=-10.1 (Mateo"1998')). We have also 
examined how the color and magnitude varies during the 
burst cycle. We find that the g — r color (which should be 
insentitive to emission lines ( West|2005 ), varies by ~ 0.2 
magnitudes during the burst and quiescent phase. The 
magnitude drops by only 0.5 magnitudes before another 
burst ignites, suggesting the galaxies should remain vis- 
ible during the quiescent phase. 

Where the models and real galaxies differ, however, is 
in the star formation rates at intermediate times. The 
models all show star formation rates that decline to the 
present day on average. However, real low mass galaxies 
tend to have nearly constant average star formation rates 
and sometimes have star formation histories that peak 
at intermediate age. This discrepancy is not surprising 
given that our models have been set up without any on- 
going source of gas accretion. Galaxies in a full cosmo- 
logical context will be capable of much more varied star 
formation histories. Our simulations suggest that very 
long separations of star formation episodes (such as in 
Carina) are unlikely to be due to feedback. Instead, the 
rec urrence of star formation must be due to the physics 
b left out of our simulations, namely interactions and on- 
going gas accretion from the cosmic web. 





0.434 



5. CONCLUSIONS 



'^•^^''Our study of isolated, spherical, collapsing low mass 
g'^g^laxies shows that star formation rates do not remain 
constant in galaxies the size of Local Group dwarfs. An 



1 Mtot is the virial mass of the halo including dark and baryoifiSective Supernova feedback mechanism is able to drive 

matter.-^ V200 is the circular velocity at the virial radius to which tgas Out of the galaxies and Create periodic episodes of 

total mass corresponds. ^ is the rms of the velocity distribution. * T^^^^^ formation. In the simple simulations presented here, 

is the stellar mass that accumulates over 10 Gyr. ° P <; p n is the time period . r j.' r n j i ■ j' 'ii i- 

corresponding to the strongest mode in a Fourier decomposition of the sM^ Star formation followed regular, periodic OSClUatlOnS. 

formation history ^ gpR Amplitude is maximum to minimum peak difffiowever, observations also show additional variabihty on 

ence of the star formation history divided by the mean star formation rate.nger time SCales including inactivity for many Gyrs. It 
b represents the median Scalo b (SFR/<SFRpast >)parameter value, j^^rd to imagine that SUCh fluctuations are pos- 

sible in the cosmological context of the Universe where 



diagrams. Star formation history recovery methods typ- 
ically use bins that increase logarithmically with time as 
shown in Figure [3] 

Overall, our simulations look remarkably similar to ac- 
tual observations. For instance, compare Figures 6 and 
7 of Young et al. ( 20071) to our 15 km s~^ star formation 



histories in I'igures [5] and [3] Like real low mass galaxies, 
the model galaxies show significant star formation at the 
present day and have not suffered from overcooling. Most 
importantly, however, they show significant variation in 
th e SFR at the present day (similar to the simulations 



jDohm-Palmer et al. 1998), GR8 (Dohm-Palmer et al. 


2002P, and Phoenix (Young et al.||2007 


i), in spite of the 



We have modeled the galaxies c olors, luminos ity and 
surface brightness using SUNRISE ( | Jonsson|2006| . SUN- 
RISE is a radiative transfer program that ray traces stel- 
lar emission through a grid that scatters light from dust, 
based on the metallicity of gas present in the simulation. 



galaxies may accrete gas sporadically over a long period 
of time and ocasionally interact with other structures in 
the cosmic web. 

Future work will discuss the outfiows from low mass 
galaxies and how star formation scales in similar galaxies 
with higher masses. 
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